Spin, Orbital and Charge Order at the Interface between Correlated Oxides 
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The collective behavior of correlated electrons in the VO2— interface layer of LaVOa/SrTiOa 
heterostructure is studied within a quarter-filled t2g-orbital Hubbard model on a square lattice. 
We argue that the ground state is ferromagnetic driven by the double exchange mechanism, and is 
orbitally and charge ordered due to a confined geometry and electron correlations. The orbital and 
charge density waves open gaps on the entire Fermi surfaces of all orbitals. The theory explains the 
observed insulating behavior of the p-type interface between LaVOs and SrTiOs. 

PACS numbers: 73.20.-r, 73.21. Cd, 75.70.Cn, 71.30.+h 




The recent progress in manufacturing and experimen- 
tal studies of heterostructures and supcrlattices based on 
the transition metal oxides has lead to the discoveries of a 
number of novel physical phenomena and new electronic 
states emerging at the interfaces [l|, 12, H Ij, la lo, LZl, la, l9| . 
A formation of a high- mobility electron gas [J] , quantum 
Hall effect [5J and, remarkably enough, even a supercon- 
ducting state y have been observed at the interfaces be- 
tween insulating oxides. At the interface between nom- 
inally non-magnetic oxides, magnetic effects have also 
been detected [7] . The physical properties of interfaces 
may largely differ from and can even be orthogonal to 
those of bulk materials, due to an "electronic reconstruc- 
tion" phenomenon [lO|. 

In correlated oxide heterostructures, electronic recon- 
struction involves not only charge, but also spin 
in particular, the orbital degrees of freedom [9|, 
since electronic orbitals are highly sensitive to the local 
environment. Such a multifaceted response of correlated 
electrons gives rise to the rich interface physics that may 
form a basis for future device applications. 

The recent work [13j added a new puzzle into this field. 
Two types of interfaces between a Mott insulator LaVOs 
and a band insulator SrTiOs have been investigated: (1) 
the V02/LaO/ri02/SrO interface with n-type polar dis- 
continuity, and (2) the LaO/y02/SrO/Ti02 with p-type 
one (formed by inserting a " metallic" SrVOa unit into the 
n-type interface). In bulk compounds, the VO2 (LaO) 
layers have a — e (-l-e) charge per unit cell, while Ti02 
and SrO layers are neutral. In such systems a polar dis- 
continuity triggers the doping of an interface layer to re- 
solve the polar catastrophe ll4|: In n-type interface the 
Ti02 layer receives a — e/2 charge while in p-type inter- 
face — e/2 charge is taken away from the VO2 layer. This 
leads to a formal valence state dP'^ of Ti and d^-^ of V at 
the n- and p-type interface layers, respectively. The resis- 
tivity measurements have shown that the n-type interface 
is metallic and p-type interface is insulating. The metal- 
lic character of Ti02 interface layer is not surprising and 
confirms existing theoretical results lid ]. However, an in- 
sulating behavior of hole doped VO2 interface layer is at 
odds with expectations and is striking, given that SrVOs 



is a good metal and already 18% Sr-doping is sufficient to 
convert bulk LaVOs into a metal, too 15]. In this Letter, 
we present a theory resolving this puzzle. In short, the 
contrasting behavior of Ti02 and VO2 interfaces origi- 
nates from their different spin and orbital structure. In- 
deed, while Ti02 layer with Ti'^+/Ti'*+ states represents 
a diluted quantum S=l/2 system (like high-Tc cuprates), 
VO2 interface layer is made of Y^+ S==l and Y'^+ S=l/2 
states - a canonical background for the double-exchange 
(DE) physics. Once spins of the VO2 layer are polarized 
by the DE mechanism, system is effectively half-filled and 
hence collective orbital and charge instabilities are trig- 
gered at the interface. We argue that these cooperative 
orderings of correlated electrons are responsible for the 
insulating character of the p-type interface. 

The Model- We describe the physics of hole doped 
VO2 layer within a multi-orbital Hubbard model for d- 
electrons [IGj on a square lattice: 



H = -J2Yl^V,^4aadjaa + U Y^ 
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The three-fold degenerate t2g states dyz, d^z, and dxy 
are labeled by orbital index a =1, 2, and 3, respec- 
tively. The Sia and n^c, — nia^ + riiai correspond to 
the spin and density of electrons in a orbital. The first 
term in H describes an electron hopping between the 
nearest neighbor (NN) sites ll7| and is diagonal in or- 
bital space. The peculiarity of a t2g system on a pla- 
nar geometry is that the orbitals dyz and d^z become 
one-dimensional (ID). They have a finite hopping ampli- 
tude only along the one particular direction, iij,i(2) = t 
for ij II x{y) and zero otherwise, see Fig. [ija). While 
the dxy orbital still forms a two-dimensional (2D) band: 
tij,3 = t. In momentum space the hopping term reads 

^s J2kaa ''^adl^^duaa, whcrc eki(2) = -2tcoskx(y) and 
Eks = — 2t(cos/i;2;4-cos fcy). The interaction part of 7J con- 
sists of an on-site intra- and inter-orbital Coulomb repul- 
sions, U and U', respectively. The latter term is further 
split by the Hund's coupling Jjj into an interorbital spin 
triplet {U' — Jh) and singlets {W + Jh), such that 5" = 1 
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FIG. 2; (Color online) Order parameters of the orbital (ti_2,3 
for dyz^xz,xy bands, respectively) and charge (S) density waves 
vs nearest-neighbor repulsion V at large U. 



FIG. 1: (Color online) (a) The hopping amplitudes of ^29- 
electrons on a square lattice, (b) The Fermi surfaces for dyz 
(horizontal lines), d^z (vertical lines) and d^y (square) elec- 
trons in the ferromagnetic state, (c) A sketch of the orbital 
and charge ordered state for U>V^t. In a symmetry broken 
state the system is insulating at any values of U and V . 



deal with interacting spinless ferniions with orbital fla- 
vors only. In this case the Hamiltonian ([T]) reduces to 



H = -V Uj,a4^d 



i,j,0! 
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state with electrons residing on different orbitals is fa- 
vored. It is this 2 Jh splitting between the high/low spin 
states that promotes a global ferromagnetic (FM) state 
by virtue of the DE mechanism in the present I/'^+/4+ 
mixed-valent system. We have also included NN repul- 
sion V which is relevant for charge ordering [18|. We 
consider a quarter-filled t2g bands {d}'^ configuration) to 
model an interface VO2 layer with formal valency V'^'^+. 
Ferromagnetism- The large on-site repulsions U 
and U' suppress the high energy charge fluctuations 
{cP,d^) ^^ (d^,d'^) and only the low energy ones 
(d^, d^) <-> {d^, d^) are allowed. In the paramagnetic and 
orbital liquid state all three bands are quarter-filled and 
there are no apparent Fermi surface (nesting) related in- 
stabilities. However, there is a potential instability to- 
wards FM ordering because the Hund's coupling favors 
the high spin state of a d^ configuration and, together 
with the hopping term, induces the DE interaction. The 
critical value of Hund's coupling at which FM instability 
appears is estimated from a vanishing determinant of a 
two component spin susceptibility for ID and 2D bands: 



1 -4JhX 

-2Jhx 1 - 2JhX 
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where x = p(0)/2 and x — p(0)/2 are uniform static 
magnetic susceptibilities for one- and two-dimensional 
bands, respectively, p(0) and p{0) are corresponding den- 
sity of states at the Fermi level. For an estimate we set 
p(0) ~ l/Wi and p(0) ~ 1/1^2, where Wi = At and 
W2 = 8t are bandwidths of ID and 2D bands, respec- 
tively. This gives a FM instability for Jh > 2{V^~ l)t ~ 
2.5i. This inequality is well satisfied for the actual pa- 
rameters Jh ^ 0.7 eV and i ~ 0.2 eV [H, [23|. We thus 



where U = U' — Jn ^ C/ — 3Jh (using well-known relation 
U' = U — 2 Jh) is an effective Hubbard repulsion in the 
FM state. Optical data in cubic vanadates [2l| suggests 
the high-spin transition at {/ ^ 2 eV ~ lOt but we will 
consider U as a free parameter and denote it below simply 
as U. Each orbital band is half-filled and the correspond- 
ing Fermi surfaces are fully nested, see Fig. [IJb). We 
discuss now the orbital and charge density waves (ODW 
and CDW) triggered by such a nesting. 

Orbital and charge density waves- Let us consider the 
instabilities towards the ODW and CDW with a modu- 
lation wave vector Q = (tt , tt). We introduce the orbital 
order parameters Tq as {uia) = n/3 -I- e'*^^'rc, where 
n — 3/2 is an average electron density. The correspond- 
ing CDW modulation is given by J^ai^ia) —n — e'*^^*(5, 
where 5 — X^q'^'q- The order parameters are calculated 
within a mean-field approach. We consider a wide range 
of U and V , in order to see how the ODW/CDW order- 
ings evolve from a weak coupling regime to the limit of 
strong interactions U ^ t. While care must be taken 
in calculating excitation spectra, the mean-field method 
gives a reliable picture of the nature of the ordered phases 
and of the zero temperature properties of interest here 
[22| . and hence is widely used to study similar problems 
(including multiorbital physics at interfaces [10|) even in 
a regime of strong correlations |23| . The ground state en- 
ergy per site is expressed in terms of order parameters: 



i^ = - J E ^k« + i ^ E ^' + J ( ^ - ^ ^ ) '5 



(4) 



where ^ka = V^L + ^l ^ith ^^ = Uto, - {U - zV)d, 
and z = 4 is a number of NNs. (The constant contribu- 
tion Eo = \Un'^ + ^zVn'^ has been dropped). Phys- 
ically, |Aq| represent band gaps. The minimization 
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FIG. 3: (Color online) Band gaps (left) and kinetic energy 
per site K (right) vs nearest-neighbor repulsion V . 
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FIG. 4; (Color online) Order parameters of the orbital (ri,2.3) 
and charge (5) density waves vs U. The dashed line marks a 
first-order phase transition (see text for details). 



of the energy E gives the coupled integral equations 
Ta = X^k^aZ-^ka solved numerically. For large U and 
V one finds T3 < T2 = — ri ~ 0.5. A sketch of a fully sat- 
urated version of the corresponding orbital and charge 
patterns is shown in Fig.[ljc). It consists of a staggered 
order of dyz / dxz orbitals, with dxy orbital being predom- 
inantly occupied at one of the sublattices. This results in 
a checkerboard charge ordered pattern of d^/d^ states. 

Shown in Fig. [2] are the ODW and CDW order pa- 
rameters as a function oi V &l U — 10 (hereafter, the 
energy unit t is used). The staggered order of dyz and 
dxz orbitals is strong and weakly affected by V (see the 
Ti and T2 curves). The strengths oi dxy orbital (ra) and 
charge {5) density waves decrease with V . Surprisingly, 
5 remains finite down to V^ = 0, i.e., a charge modula- 
tion is present even in the case of local interactions only. 
A similar effect but with different mechanism has been 



found within the two-orbital model for manganites [24 1. 
Physics behind this unusual picture here is as follows. In 
the limit U ^ t the orbital order parameters ri^2 for ID 
bands are nearly saturated and can be expanded in pow- 
ers of t/U. The energy (|4]) is then expressed in terms of 
rs only: E = -\U - ^ + 2 Jt| - i Ek V^L+WW, 
where J = At^ /U . The minimization of E gives a finite 
T3 and non-zero CDW order parameter 5 ~ T3 — 0{J/U), 
because of a singular response of the nested dxy Fermi 
surface at half-filling. 

The orbital and charge density waves with Q = (tt , tt) 
modulation wave vector induce the gaps on the entire 
Fermi surfaces of all three bands and drive the system 
into the insulating state. In Fig. [3] the band gaps are 
plotted as function oi V at U = 10. At 1/ = 0, the gaps 
in the dyz and dxz bands (|Ai| and IA2I, respectively) 
are large since they scale as U in the limit U ^ t. While 
the gap of dxy band is controlled by an effective coupling 
constant oc J. In the "weak-coupling" limit J ^ t {i.e. 
strong-coupling U ^ t limit in conventional language) it 
is exponentially small, jAaj ~ texp[— i/4J] at V^ = 0. At 
large V, the expected [Aaj ^ 2V scafing is observed. 

To complete our analyses, in Fig. 4 the dependence 
of the order parameters on U is presented for T^ = 0. 



The dashed line marks a first-order phase transition at 
around Uc — 2.4 from the phasc-I, sketched in Fig. [Ijc), 
to the phase-II. In the latter dyz and dxz orbitals pre- 
dominantly occupy one sublattice while the density of 
dxy orbital is higher at the other one. The ODW order 
parameters for dyz and dxz bands (ri and T2) are con- 
trolled by U and monotonically decrease with U. The 
non-monotonic behavior of dxy orbital and CDW order 
parameters is explained as follows. The effective coupling 
constant controlling them vanishes in the limits of small 
as well as of large U: being of the order of U for U <^ t 
and ex t^/C/ for U ^ t. For realistic values of model pa- 
rameters the phase-I, sketched in Fig. [TJc), is the ground 
state of the system. We emphasize that the system is 
insulating for any finite values of U and V. 

Thus, we arrived at rather unusual situation where DE 
driven FM and insulating states coexist 25| and, more- 
over, are closely interrelated. In fact, the ODW and 
CDW states are stable only if FM correlation length is 
large. In the DE system, a kinetic energy of electrons 
defines a stiffness of FM order. Along the same line we 
estimate a FM coupling of neighboring s = 1/2 and S = 1 
spins in the charge ordered state: Jfm — ^/[2S'(2s + 1)], 
where K is a kinetic energy per site. Fig. [3] shows K 
as a function of V. Considering a moderate value of 
y = 2i ~ 0.4 eV, we find Jfm ^ 40 meV. This sug- 
gests an onset of FM correlations at fairly high temper- 
atures and explains insulating behavior of VO2 interface 
in the experiment [13[. Further, we predict a transition 
to a metallic state and large magnetoresistivity effects at 
higher temperature when FM correlations are reduced. 
Apart from transport measurements, magnetic x-ray and 
optical studies may provide a crucial test for the theory. 

The present work motivates an interesting idea of a 
superlattice depicted in Fig. [5l Here, insertion of SrO 
planes into LaVOa can be viewed as a "spatially corre- 
lated doping" that generates two V'^-^+02 planes - each 
midway between Sr^+0 and La'^+O layers - forming a 
ferromagnetic bilayer. Spins of different bilayers weakly 
couple antiferromagnetically (ferromagnetically) if their 
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FIG. 5: (a) V^'^^02 bilayer formed by replacement of a 
LaO (001) layer of LaVOs by SrO. (b) Spin structure of the 
yS.s+Q^ bilayer. Large (small) circles and spins indicate cites 
where V '^dxy-jdyz^ (V '^d^x^) configuration is favored. The 
layers are coupled ferromagnetically via the DE between S=l 
and S=l/2 states, (c) Periodic sequence of V'^'^'^02 bilayers. 
Spin coupling through the intermediate V^"*'5' — 1 ions of 
VO2 layer is antiferromagnetic (as shown here) if the charge 
ordering patterns of different bilayers are in-phase. 



CDW-ordering patterns are in-phase (out-of-phase) [26|. 
Such a direct link between charge and spin structures sug- 
gests a magnetic control of charge sector and vice versa, 
e.g., a relatively weak magnetic field may lead to the 
CDW phase-shift. In fact, the proposed superlattice is 
similar to the bilayer ruthenates [27| and manganites [28| , 
which have the same magnetic structure as in Fig. EJc) 
and show large magnetoresistivity and spin- valve effects. 
We notice also that the bilayer coupling may oscillate in 
sign as the number of intermediate VO2 planes is varied, 
provided their V'^"'" spins stagger along the c axis as in 
YVO3 at low-temperature [21 1. 

To conclude, we have studied a Hubbard model for a 
quarter-filled ^23 bands on a square lattice. Due to a con- 
fined geometry at the interface, at large enough Hund's 
couphng (Jh ^ 2.5i) the correlated electrons develop a 
peculiar insulating ferromagnetic ground state accompa- 
nied by the orbital and charge density waves. This pro- 
vides a natural explanation for an insulating behavior 
of the p-type LaVOa/SrTiOa interface. The experimen- 
tal and theoretical studies of superlattices like in Fig. [5] 
where a complex electronic reconstruction takes place co- 
herently over many interfaces remains a future challenge. 
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